Electrical Breakdown in a V2O3 device at the Insulator to Metal Transition 
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We have measured the electrical properties of a V2O3 thin film micro bridge at the insulator metal 
transition (IMT). Discontinuous jumps to lower voltages in the current voltage characteristic (IV) 
followed by an approximately constant voltage progression for high currents indicate an electrical 
breakdown of the device. In addition, the IV curve shows hysteresis and a training effect, i.e. the 
subsequent IV loops are different from the first IV loop after thermal cycling. Low temperature 
scanning electron microscopy (LTSEM) reveals that the electrical breakdown over the whole device 
is caused by the formation of electro-thermal domains (ETDs), i.e. the current and temperature 
redistribution in the device. On the contrary, at the nanoscale, the electrical breakdown causes the 
IMT of individual domains. In a numerical model we considered these domains as a network of 
resistors and we were able to reproduce the electro-thermal breakdown as well as the hysteresis and 
the training effect in the IVs. 



MAIN PAPER 

Stoichiometric V2O3 is a strongly correlated material 
that undergoes a first order insulator to metal phase 
transition (IMT) from an antiferromagnetic insulating to 
a paramagnetic metallic phase. Because the IMT causes 
a change in the resistivity of several orders of magnitude 
novel devices based on the IMT and their applications 
are actively investigated [1 and it is in focus of current 
research whether there is a voltage driven IMT in 
strongly correlated materials [2-6 . If this hypothesis 
holds, a micro bridge fabricated of a strongly correlated 
material would experience a dielectric breakdown for 
high bias voltages. However, because of the large 
resistivity change, self-heating in such a device can be 
strong enough to cause an electro-thermal breakdown, 
i.e. an current and temperature redistribution [7[[8]. In 
order to investigated the electrical breakdown of a V2O3 
micro bridge, we used a very unique low temperature 
scanning electron microscope (LTSEM) to image the 
metallic and insulating phases and developed a numer- 
ical model to simulate the electrical device properties. 
Here we show that the electrical breakdown over the 
whole device is due to the formation of electro-thermal 
domains (ETDs), and that the distribution in the IMT 
temperature as well as the hysteresis in the RT of 
domains at the nanoscale are significantly infiuencing 
the current voltage characteristic (IV). 

The V2O3 film was grown by rf-sputtering on a r-cut 
sapphire substrate at a temperature of approximately 
750 °C . XRD revealed that the polycrystalline film grows 
textured under these conditions. Optical contact lithog- 
raphy was used for patterning. A 200 jiia wide V2O3 
bridge was etched using reactive ion etching and gold 



electrodes were evaporated on top, refer to Fig. [T] (a) for 
an SEM image of the device. 

The same experimental set up was used for LTSEM 
imaging and IV characterization. The LTSEM is a con- 
ventional, state of the art SEM equipped with a LN2- 
cryostat. The sample is mounted in a vacuum on a cold 
plate and the sample temperature can be varied from 
room temperature to approximately 80 K. Due to the 
large working distance required by the IR-shielding of 
the cryostat the smallest spot size of the electron beam 
is about 300 nm. A Keithley 2400 current source and 
a conventional preamplifier with a National Instruments 
measurement card were used for the two point electrical 
measurements. 

To acquire a LTSEM image, the sample is scanned with 
a blanked electron beam and the electric response /SV of 
a simultaneously applied four probe or two probe mea- 
surement is mapped in a so called voltage image of the 
sample [9l[T0]. 

We have measured the current-voltage characteristics 
of a 200 /im wide and 10 /im long V2O3 bridge (Fig. [l^) 
for the base temperatures T5 = 145 K and T5 = 155 K, 
respectively. Due to the IMT there is a resistance change 
by more than three orders of magnitude in a 20 K tem- 
perature interval (Fig. ^p) and a pronounced hysteresis, 
i.e. the heating curve is shifted to higher temperatures 
by approximately 5 K. 

To ensure that the V2O3 film was completely in the in- 
sulating state the device temperature was reduced to ap- 
proximately 80 K and then heated to the intended base 
temperature before electrical measurements were started. 
In the following, we refer to this procedure as "thermal 
cycling" . 

Fig. Ic shows an IV of the pristine device (after ther- 
mal cycling) acquired at a base temperature T5=145K. 
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FIG. 1. (Color online) Electrical breakdown, (a) SEM image 
of the device. The area under investigation is indicated by a 
red rectangle. A 200 /xm wide V2O3 stripe runs vertically 
indicated by two white dashed lines. The two gold electrodes 
at the top form a 10 /am wide gap. The current flows between 
those two electrodes vertically, (b) Temperature dependence 
of the device resistance. Current voltage characteristics (IVs) 
at two different base temperatures. The black arrows indicate 
the sweep direction. The overall progression of the electrical 
breakdown is indicated by red dashed arrows, (c) and (e) are 
IVs of the pristine device after cooling to 80 K. (d) and (f) 
are subsequent IVs to (c) and (e), respectively. 



Starting at the origin of the graph, the IV shows an al- 
most linear dependence with a slight upwards bent. At 
a current value of 0.6 mA there is a jump from 10 V to 
7V followed by several small jumps and a big jump to 
3.8 V. The curve continues almost vertically with a pos- 
itive slope and several sawtooth like discontinuities of 
different sizes. The down sweep curve is different from 
the up sweep, i.e. there is a hysteresis. Immediately after 
the first IV curve, without changing the base tempera- 
ture, a second IV was acquired (Fig. [T]F). Note, that the 
IV of the second sweep is very similar to the first one ex- 
cept that the maximum voltage is significantly reduced. 
When the device is thermally cycled, the old state (first 
sweep) is restored. Hence, there is a training effect. 
The curve acquired at the first sweep at 155 K (Fig. [l^) is 
different from the one at 145 K. Starting at the origin, the 
IV is almost linear. After 1mA and 3.5 V it progresses 
with a negative slope and saw tooth like discontinuities 
of different sizes. The down sweep curve is rounded with 
some very small discontinuities and the hysteresis is much 



more pronounced than at 145 K. In the second IV (Fig{T]F) 
the down sweep is similar to the down sweep of the first 
IV, but the up sweep curve is round and there are only a 
few small discontinuities. A comparison between the IV 
of the first and second sweep reveals, that the training 
effect at 155 K leads to a qualitative change of the up 
sweep curve. The consecutive IVs do not change signifi- 
cantly after the second sweep at both base temperatures. 
Moreover, we have observed that the details of the IV de- 
pend on the sweep velocity. If the current is swept faster, 
there are more and smaller jumps, but the overall shape 
of the IV stays the same. 

At T5=145 K a big jump in the IV followed by a vertical 
progression indicates an electrical breakdown of the de- 
vice. The diagonal progression with a negative slope and 
small jumps in the first IV at T5=155K implies that the 
electrical breakdown now evolves via stable intermediate 
states. From the hysteresis in the IVs and the training 
effect it can be inferred that there is a memory effect in 
the V2O3 device, i.e. the device properties depend on 
its history. Similar effects were reported in VO2 devices 
PJJ . The slight influence of the sweep velocity on the IV 
might be caused by a slowed down relaxation of a small 
fraction of the film due to the spread in the IMT tem- 
peratures of individual domains [12 . 
The LTSEM images at a base temperature of 145 K, are 
discussed below [Fig. |2] . In order to relate the images to 
the current voltage characteristic, the IV of the second 
sweep in Fig. [l] is depicted. We have discovered that 
the imaging process changes somewhat the resistance of 
the device. In Fig. |2] (a) the bias points before and af- 
ter imaging are indicated by red arrows. In the images 
at bias points below the first big jump (b-e) a cluster of 
dots appears in an approximately 10 /im wide section of 
the device indicated by a dashed yellow rectangle in im- 
age (f). In the images at bias points above the first jump 
(f-k) a bright stripe appears whose width increases with 
increasing currents. 

A series of LTSEM images acquired at a base tempera- 
ture of 155 K are shown in Fig. [3) If the bias is increased, 
submicron sized dots of different brightness appear. The 
majority of them is clustered in the same small section 
of the device, where at T5 = 145 K the electrical break- 
down occurs. For increasing currents the density of these 
dots increases and they merge. When an imaging scan is 
repeated at the same bias current, the images look more 
or less the same and additional bright spots appear at 
different positions. 

In the supplement, the LTSEM response AV to the mod- 
ulated electron beam is estimated. This shows that A V is 
proportional to the temperature derivative of the conduc- 
tivity close to the electron probe Because ^ in the 
insulating phase is several orders of magnitude smaller 
than in the metallic phase, a large response AV is pre- 
dominately caused by the metallic domains. The bright 
sub microns sized spots, which saturate the signal in the 
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FIG. 2. (Color online) LTSEM volt age- images acquired at 
different bias currents at a base temperature of 145 K. (a) IV 
of figure 1 (c) . The start points and the end points before and 
after imaging are indicated by red arrows, (f) image of the 
entire device, (b^ — k^) images with a reduced field of view 
(indicated by a yellow dashed rectangle in f). 



LTSEM images, might be due to current redistribution. 
In the LTSEM image series at T5=145K (Fig. [2| a 
metallic filament appears (electro-thermal domains) like 
reported in VO2 devices [7l[8]. In the T5=155K series 
no filament appears, but the device becomes gradually 
metallic in a small section of the device when the current 
is increased. 

In order to investigate how the IMT transition of 
domains at the nanoscale are influencing the electrical 
breakdown, we have developed a numerical model (simi- 
lar to the one used in [13 ), in which we represent domains 
with different IMT temperatures by a 20 x 400 resistor 
network (see supplement). For each domain the same 
hysteresis ATc in the RT-dependence was assumed, but 
for the IMT temperatures Tc a Gaussian distribution was 
used (Fig. [4] a). The hysteresis as well as the median and 
the standard deviation of this Gaussian distribution were 
obtained by optimizing the simulated RT-curve with re- 
spect to the measured RT (Fig. [4]b). Then these values 
(ATc = 8K, Tc = 163 K and RMSD^ 3.16) were used 
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FIG. 3. (Color online) LTSEM voltage images acquired at 
different bias currents at a base temperature of 155 K. (a) 
IVs of figure 1 (e-f). The start points and the end points 
before and after imaging are indicated by red arrows, (b, e, 
and k) images of the entire device, (b^^ — k^) images with a 
reduced field of view (indicated by a yellow dashed rectangle 
in b,e, and k). 



to simulated the IVs as well as the temperature, cur- 
rent and voltage distributions for the base temperatures 
n = 145 K (Fig. |4]c-e) and n = 155 K (Fig. |4]f-h). The 
simulated current and voltage distributions are shown in 
the supplement. 

The hysteresis and the training effect in the IVs could 
be reproduced. The simulation at T5 = 145 K shows the 
formation electro-thermal domains like in the LTSEM- 
images (Fig. [2|. There are significant differences be- 
tween the simulation and measurements at T5 = 155 K: 
The simulation clearly shows the formation of electro- 
thermal domains (filaments) , while the IV progresses ver- 
tically (Fig. |4]f-g). In the LTSLM images no filament 
appears and the IV progresses diagonally (Fig. [3|. 
Conclusion: The hysteresis and the training effect in the 
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FIG. 4. (Color Online) Numerical model, (a) Assumed resistance vs. temperature dependence for an individual domain. Tc 
denotes the IMT temperature (heating curve). ATc is the temperature difference by which the cooling curve is shifted, (b) 
Simulated resistance vs. temperature dependence. The crosses are simulated points and the continuous line is the measured 
curve, (c) Simulated IV of the pristine device (base temperature Tb = 145 K). (d) Consecutive IV. (e) Temperature distributions 
in the device at different bias points (marked in (c) and (d) with Greek letters), (f) Simulated IV of the pristine device (base 
temperature Tb = 155 K). (g) Consecutive IV. (h) Temperature distributions in the device at different bias points (marked in 
(f) and (g) with Greek letters). 



IV are the results of the RT hysteresis as well as the 
distribution of IMT temperatures of nanoscaled domains 
in the polycrystalline thin film. It is plausible that the 
RT hysteresis originates from the structural bistability 
reported earlier p^HTB] and the spread in the transition 
temperatures is most likely caused by differences in the 
strain in the polycrystalline film [121 III] • 
At a base temperature of 145 K the overall IV and the 
ETDs (filaments) depicted in Fig. [5] (f-k) were repro- 
duced in the simulation, clearly indicating an electro- 
thermal breakdown of the device. 

On the contrary, at a base temperature of 155 K the 
electrical breakdown, with the small metallic domains 
clustered in a small section of the device and the diag- 
onal progression of the IV (Fig. [3|, is very atypical for 
an electro-thermal breakdown. Generally, in a system, 
where the electro-thermal bistability is caused by a large 
decrease in the device resistance with increasing temper- 
ature the electro-thermal domain walls are parallel to the 
current direction, i.e. in the sample under investigation, 
electro-thermal domains have always the shape of a fila- 
ment. After electro-thermal domains have nucleated, the 
IV progresses vertically, while the hot electro-thermal do- 
mains increase in size j8l|T9] (see also the supplement). 
We provide three different explanations for the electri- 
cal breakdown via stable intermediate states at a base 
temperature of 155 K. First, it is possible to stabilize an 
electro-thermal bistable device in intermediate states by 
using a load resistor, which forces the IV to progress 



along a diagonal load line during the electro-thermal 
breakdown [20j. The device under test has a relatively 
high contact resistance of approximately 80 Q. If one as- 
sumes that this contact resistance is nonohmic and has a 
tunneling characteristic, the sawtooth like discontinuities 
and the diagonal progression of the IV at 155 K can be 
explained. Second, an intrinsic shunting of small domains 
above a certain threshold voltage that could be caused by 
Landau-Zener tunneling [2]-[5] might have a similar effect. 
Third, in a theoretical model considering voltage induced 
switching two different electrical breakdown mechanism 
were predicted - bolt like and percolative switching [21]. 
According to this results one might interpret the electri- 
cal breakdown at a base temperature of 145 K as bolt 
like and at 155 K as percolative switching. In this model 
thermal coupling within the thin film and to the sub- 
strate was not considered. When we explicitly included 
voltage induced switching in our model, we were never 
able to observe percolative switching, but we always ob- 
served the formation of a filament. The only effect was 
to shift the thermal breakdown to lower voltages. This is 
quite plausible, because the sample heating effects occur 
at lower voltages, in this case. From the simulated volt- 
age distribution in the supplement a threshold voltage for 
a dielectric breakdown below 20kV/cm can be inferred. 
This value would be very small. Therefore, we consider 
the first explanation as the most likely scenario. 
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SUPPLEMENT 

Self-heating, electrothermal bistability and the 
formation of electrothermal domains (ETDs) 

Self-heating, thermal bistability and the formation of 
electrothermal domains (ETDs), are discussed in the fol- 
lowing section. The primary cause for these phenom- 
ena is a strong change in the resistivity vs. temperature 
characteristic of the material used for the device under 
investigation. These phenomena can be observed in large 
number of different kind of devices: There is the hot-spot 
in a superconducting microbridge [191 1221 [22] , the pinch 
caused by thermal breakdown in a negative temperature 
coefficient thermistor (NTC) [211 ES], the domains in a 
two valley semi-conductor used in a Gunn-diode [18] , the 
hot spot in a large BSCCO-mesa used as a THz-emitter 
[26H29] or the filaments in a VO2 device [T] |30] at the 
metal-insulator transition, to name only a few examples. 
In the following, only the basics of this topic necessary 
to understand the supplemented paper are provided. A 
more elaborate discussion can be found in the review ar- 
ticles [nuin]. 



Self-heating 

In a device with finite electrical resistance the Joule 
heating leads to a rise in the device temperature. This 
temperature increase depends on the thermal coupling 
of the device to the environment. If the device resis- 
tance is temperature dependent (as it is usually the case), 
this self-heating effect causes a deviation in the current 
voltage characteristic (IV) from the linear (Ohmic) pro- 
gression. In other words, there will be bending in the 
IV. If the device resistance changes strongly with tem- 
perature, this bending can be very pronounced as it is 
demonstrated in the following. 

If the temperature dependence of the device resistance 
and the thermal coupling are known, and a uniform tem- 
perature distribution is assumed the IV can be parame- 
terized using a method similar to the graphical method 
described in [31. 

The excess temperature 1^ of the device is determined by 
Ohm's law and Newton's law of cooling: 



V = R{i9)I; VI = A{i3), 



(1) 



where A{d) is the heat transfer coefficient. The parame- 
terization of the IV is 



V = ^JR{'^) A{;d) and / : 



wr 



(2) 



This method is now applied to the V2O3 device discussed 
in the article. On the assumption that the thermal con- 
ductance of the device is solely determined by the sap- 
phire substrate and "edge cooling" can be neglected the 
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FIG. 5. Influence of self- heating on the current voltage char- 
acteristic of the V2O3 device: (a) Temperature dependence of 
the device resistance, (b) Thermal conductivity of sapphire 
vs. temperature, (c) and (d) calculated current voltage char- 
acteristics of the device, assuming a homogenous tempera- 
ture distribution. The base temperatures are Th — 145 K and 
Th = 155 iT, respectively. Inset (e) heat transfer coefficient 
via the sapphire substrate vs. excess temperature assuming 
a base temperature Th — 145 K. Inset (f ) IV of figure (d) at 
a larger scale. 



heat transfer coefficient can be calculated using the ther- 
mal conductivity n of sapphire 



LW 



n(T) dT 



(3) 



Here L is the length, W the width and H the height of 
the device and T5 is the base temperature. 

Figure [s] c) and d) show the parameterized IV at the 
base temperature T5 = 145 K and T5 = 155 K, respec- 
tively. H = 0.5 mm, W = 200 /im and L = 20 /im were 
used. For current values below 1 mA the calculated IVs 
of this supplement are similar to the measured IVs pre- 
sented in the main paper. The IV with a large current 
scale (inset [Fig. jsjf)]) has a very pronounced S-shaped 
characteristic. 



Electrothermal Bistability 

In the preceding section it was demonstrated that in 
a device, with a strong resistance decrease in a small 
temperature interval, self-heating leads to an S-shaped 
IV curve, cf. [Fig. jsjf)]. If a voltage source is used to 
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FIG. 6. Formation of electrothermal domains, according to 
[24j : (a) NTC clamped between two metal pins, (b) Thermal 
equivalent circuit: Ai and A2 denote the thermal coupling to 
the heat sink and A12 the thermal cross coupling, respectively, 
(c) and (d) IV curve corresponding to the area 1 and area 2 
of the device, (e) IV of the whole device. 



bias such a device, then the system is (electrothermally) 
bistable, i.e. there exists a bias interval in which the 
system can rest in two different bias states: The one 
state has a low bias current and a low device temperature 
(denoted in [6] (c) and (d) with /3i and ^2)^ while the 
other state has a high bias current and a high device 
temperature (denoted in[6](c) and (d) with 71 and 72). 
Note, that for a device with a strong increase in resistivity 
the self-heating is causing an N-shaped IV and voltage 
and currents are changing their roles, but this supplement 
is restricted on the case with an S-shaped IV. 



Formation of Electrothermal Domains 

Electrothermal domains can form in a device with weak 
thermal coupling between different parts. According to 
[32] one of the first publications in the field of ETDs 
was the work of E. Spenke on the thermal breakdown 
of negative temperature coefficient thermistors (NTCs) 
p4( [25] . He discussed a NTC disk fixed between two 
metal pins, which are in contact with a thermal reser- 
voir at a constant temperature, [Fig. [6]^ a)]. Note, that 
the resistance of the pins and the thermal reservoir is 
considerably smaller than the resistance of the disk. In 
a "Gedankenexperiment" , E. Spenke divided the NTC- 
disk into two domains, or in other words he considered 
the disk consisting of two NTCs in parallel, which are 
thermally coupled via the metallic pins. He inferred that 
the system (under current bias) can stabilize itself in the 
bias interval of negative differential resistance, if one of 
those NTCs switches to a bias point in the IV that corre- 
sponds to a high bias current and high temperature (for 
instance bias point 71) and the other NTC switches to a 
bias point with low current and low temperature (for in- 
stance bias point ^2)- This is equivalent to the formation 
of two areas with different current density and temper- 
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FIG. 7. a) Scheme of the device under investigation (top 
view), b) sketch of the specific conductivity vs. temperature. 



ature, which are called electrothermal domains (ETDs). 
Note, only if the ETD walls are parallel to the current di- 
rection, the ETDs can be considered as parallel resistors 
and the argumentation of the "Gedankenexperiment" can 
be applied. 

The formation of the electrothermal domains is causing 
a kink or a discontinuity in the IV. If the S-shape in 
the IV of a device with uniform temperature distribution 
is very pronounced than this discontinuity is very pro- 
nounced, too. After the formation of the ETDs the IV 
progresses almost vertically indicating that an increase of 
the current is increasing the size of the "hot" electrother- 
mal domain with a high current density while the voltage 
over the device stays constant. 



Low temperature scanning electron microscopy 
(LTSEM) 

In the following we derive a relation between the LT- 
SEM voltage signal AV(xo, ^0) and the local conductance 
per area g of the V2O3 sample. In LTSEM a periodically 
blanked focused electron beam (using /5 = 13.3 kHz, 
acceleration voltage 10 kV, beam current 100 pA), is 
scanned across the sample surface in the (x, ?/)-plane. 
When the electron beam is positioned at coordinates 
(xo,^o), it causes local heating and thus local changes 
in temperature-dependent parameters such as the local 
conductivity of the junction. This can be detected by 
means of lock-in technique as a beam-induced change 
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AV{xo,yo) of the voltage V across the junction, which 
is biased at a constant current /. The beam current also 
adds to the bias current density in the beam spot around 



{xq^Vo)' For measurements reported here, the beam cur- 
rent density is several orders of magnitude smaller than 
the typical transport current densities. Therefore, this 
effect will be ignored here. 



If the electron beam is on, the local temperature increases by 5T resulting in a temperature distribution T{x — 
xo^y — yo) = To -\- ST{x — xo^y — yo)^ where Tq is the local temperature of the undisturbed junction. The temperature 
profile created by the beam typically has Gaussian shape 



6T{x - xo,y - yo) ATexp < - 



and AT ?^ 0.1 - 1 K and a ^ 0.3 - 2 /im. 



{x - xo)^ + - yoY' 



2^2 



(4) 



The geometry of the V2O3 device is a microbridge of 
length L along x, with a width W along y and a thick- 
ness D along z [Fig. [7|(a)]. Along z the incident electron 
beam heats the bridge almost homogeneously and we as- 
sume that the sample specific conductivity g{x^y) does 
not depend on z. The bias current is applied in x direc- 
tion. Further, since W ^ L we assume that the currents 
inside the bridge strictly flow in x direction, i. e. we 
neglect voltage drops in y direction and current redis- 
tribution effects. With an electron beam illuminating a 
spot at the position (xo,yo) the beam induced voltage 
change is given by 



AV{xo,yo) = H^on - ^off), 



(5) 



where / is the total device current and i^on/off is the 
resistance when the laser beam is on or off, respectively. 
We divide the device in N slabs oriented in the (y,z) plane 
with thickness Ax. The i-th slap with the coordinate Xi 
has the conductance Gi and the electron beam causes a 
small change in conductance SGi. The device resistance 
for the electron beam off state can be written as i?off = 



N 



= 1 Gi 



and using a Taylor expansion 
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(6) 



From [5] and [g] follows 
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5G^ 
G? 



(7) 



We approximate the conductance of the i-th slab by av- 
eraging over the device 



Ax i^off 



SGj can be written as 



6G, 
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5g{xi,y,xo,yo)dy, 



(8) 



(9) 



where Sg{x^y^Xo^yo) is the local change of conductivity 
cause by an electron beam illuminating a spot at the 
position {xo.yo)- 



We obtain for the beam induced voltage change 
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y,xo,yo) Axdy = 



ID 
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/ 6g{x,y,xo,yo)d{x,y) (10) 

J A 



The quantity to be evaluated further is 
J^dg{x^y) d{x,y). In analogy with the assumed 
R{T) dependence of an individual domain, see inset 
in [Fig. 5] of the main paper, we consider g{x^ y) of 
the domain as a function of T, which is characterized 
by three regions, see Fig. [t] (b)]. At temperatures Tq 
far away from the insulation-to-metal transition (IMT) 
g{x,y,To) is a unique function with g^ for the insulating 
state and gm for the metallic state, respectively. At the 



IMT the curve is hysteretic. 

The e-beam induced change in conductance of an 
individual domain is qualitatively discussed in the 
following: Tq is the temperature of the domain when the 
e-beam is off and ST is the temperature increase of the 
domain when the e-beam is on. According to [Fig. [t], 
if To > Ti (where T^ is the insulator to metal transition 
temperature) the change of the domain conductance 
is proportional to the slope If Tq < Tj and 
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To -\- 5T < Ti^ the change of the domain conductance 
is proportional to the slope ^ at Tq. If Tq < Ti and 
< To + (5T, the illumination with the electron beam 
causes an irreversible switching of the domain. This 
mechanism explains why the LTSEM imaging process is 
influencing the bias over the device as discussed in the 
main paper. After the domain has switched the beam 
induced change of the conductance is proportional to 
Finally, if Tq < (where is the metal to 
insulator transition temperature) and < Tq + 5T^ the 



domain switches reversible between the insulating and 
the metallic state and the beam induced change of the 
domain conductance is proportional to Qm — Qi ^ gm- 
Due to the lock-in technique this response can be 
detected in LTSEM in contrast to the irreversible 
switching. To summarize, there are three different 
response mechanism of an individual domain: Metal- 
lic, insulating and the reversible switching, which is 
denoted with the index h (for hysteresis) in the following. 



Because the device temperature locally varies and different domains have different IMT temperatures, the device is 
divided in three areas A = Am -\-Ai-\-Ah. Note, that the areas are not necessarily connected. For instance the device 
could consist of small metallic "puddles" embedded in an insulating matrix. 
We rewrite Sg{x, y) d{x^ y) as 



/ 5g{x,y)d{x,y) = 

J A J A, 

r [dgi{x,y)^ 

JaA d: 



|to ^T{x -xo,y- yo) 



dgmjx.y) 
dT 



|to ST{x -xo,y- yo) 



d{x,y) 



(11) 



d{x,y)^ / [gm{x,y) - gi{x,y)]d{x,y). 

J Ah 



The first term on the right hand side represents the re- 
gions of the sample which remains in the metallic {gm) 
state and the second term represents domains that re- 
main in the insulating state {gi). The third term repre- 
sents the domains, which reversibly switch between the 
insulating and the metallic state upon illumination. We 
may assume that TV domains in a radius ~ cr switch be- 
tween the insulating and the metallic state. Let the area 
of the domain n be An- Then the third term reduces 
to Y^n^n [gm{Tm) - ^i(^i)] 6( |ro - Tn | /cr) , where ro is a 
vector to point (xq, ^o) and points to the center of the 
domain n. 6(^) shall be 1 if ^ < 1 and otherwise. If 
is on the order of or smaller than A^ the domains which 
switch will be clearly distinguishable and A/" = or 1. For 



cr^ An and not well separated switching domains one 
will see some blurred signal where the switching domains 
cannot be resolved clearly. 

In the following we calculate and discuss the expected 
different responses of Eq. (11) for a device at 145 K 



at a bias 10 mA and 3.5 V which corresponds to image 
k in Fig. [3] of the main paper. For the insulating 
phase gi and ^ can be established by means of the in- 
set in Fig. [5] of the main paper. At T = 145 K we get 



gi 



and ^ 



1.37 



Because the device 



we discussed in the main paper has a high contact resis- 
tance, we used a four-terminal measurement on a device 
with comparable V2O3 thin film properties to determine 
gm - 1.43 • 10^ and ^ ^ -310, ' 



For the beam induced voltage change due to switching of domains we finally find 



AVh{xo,yo) ^ -IRl^^yAngm(To)Q{\ro - r^l/a) 



n 

-0.05 V, 



(12) 
(13) 



with N =1 and A^ (0.5 x 0.5) jum^ 



r 



This expected value exceeds our measured signal by sev- 
eral orders of magnitude and is even higher for lower 
bias currents. According to the simulations of tempera- 
ture distribution [Fig. 6-7 of the main paper], the local 
temperature within an ETD (Tetd ^ 210 K) is clearly 



above T^. Furthermore, since the electron beam gener- 
ates a local increase in temperature of AT = 0.25 K and 
the hysteresis has an amplitude of T^ — T^ = 8 K AT, 
we rule out a periodic switching of domains. 
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For the induced voltage change AVi^m due to domains that remain in the insulating state or metallic state respec- 
tively, we find by using Eq. Q: 

AV,,^{xo,yo) « -IRl^^AT l^^{xo,y,)^ (14) 

where ^ ^^y^ (xq, yo)^ is the convolution of the local value of ^^y^ (xq, y^) with the electron beam induced temperature 
profile. 



In order to obtain an upper estimate for the induced volt- homogeneously heated up by the electron beam and that 
age change AV, we assume that an area Ag = ^(Ax)^ is the entire area Ag is either metallic or insulating. 



J 



In this case equation 14 reduces to 



Ayi,™(xo,yo) « -/i?^ff^ATA,^(xo,yo) (15) 



By considering the resolution of the LTSEM images in the 
main paper we can estimate Ax ^0.5 /im. We assume a 
beam induced temperature increase AT ^ 0.25 K. This 
value is reasonable. For instance in [33] AT = 0.2 — 0.4 K 
was obtained by considering a totally different LTSEM 
response mechanism. Finally, we find /SVi ^ — 82 nV for 
the insulating phase (which may be neglected due to its 
smaller magnitude) and /SVm ~ 19 /iV for the metallic 
phase. From this we infer that mainly the metallic phase 
causes the response AV in the LTSEM images. 
Note, that the preceding estimate is an upper limit. The 
signal is reduced, if not all domains within the electron 
beam induced temperature distribution are in the metal- 
lic state. This explains the variation of brightness in 
the LTSEM images. In addition, current redistribution 
effects are neglected in the preceding model. But these 
current redistribution effects can be considerable in a net- 
work of fine conducting metallic filaments embedded in 
an insulating matrix. This might be the reason for the 
very bright spots in the LTSEM images of Fig. 4 in the 
main paper. 



Computer Model and Additional Simulation Data 

In order to simulate the Joule heating in the device, we 
have developed a computer model, in which the domains 
are represented by a diamond shaped resistor network 
[Fig. [l2|a)] (like in [21]) consisting of M x diamonds. 
Note, that the program we used is similar to the one used 
in [13]. For the calculations, the electrical properties as 
well as the temperature of these resistors are represented 
by 4 X M X dimensional vectors, and for the visualiza- 
tion, we use two dimensional matrices [Fig. [T2|b)]. To 
simulate the IV, the program starts with an input cur- 
rent /in = and a homogenous temperature distribution 
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FIG. 8. Simulation of the V2O3 device with a base temper- 
ature of 145 K. (a) IV of the pristine device, (b) consecutive 
IV. (c-k) current distribution in the device at different bias 
points (marked in Fig. a and b). 



T equals to the base temperature T5. For every resistor a 
temperature dependence of the resistance like in the in- 
set of Fig. 4 in the main paper with individual transition 
temperatures is used. The conductance vector G is cal- 
culated using the temperature distribution T. We have 
developed an algorithm that creates a set of 4 x M x 
independent Kirchhoff equations corresponding to the re- 




FIG. 9. Simulation of the V2O3 device with a base temper- 
ature of 145 K. (a) IV of the pristine device, (b) consecutive 
IV. (c-k) voltage distribution in the device at different bias 
points (marked in Fig. a and b). 



FIG. 11. Simulation of the V2O3 device with a base temper- 
ature of 155 K. (a) IV of the pristine device, (b) consecutive 
IV. (c-k) voltage distribution in the device at different bias 
points (marked in Fig. a and b) 
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FIG. 10. Simulation of the V2O3 device with a base temper- 
ature of 155 K. (a) IV of the pristine device, (b) consecutive 
IV. (c-k) current distribution in the device at different bias 
points (marked in Fig. a and b) 



sistor network 

KM(G)F = 4. (16) 

We solve for the voltage distribution 

V = KM-14, (17) 

and calculate the current and power distributions 

T=G-V, P = I-V. (18) 

With the power distribution the temperature distribution 
can be recalculated. For the thermal cross coupling (cou- 
pling parameter Ac), we consider only coupling between 
nearest neighbors. For instance, for a network consist- 
ing of only 9 resistors (depicted in Fig. [l2|c)) we get a 
thermal coupling equation 

Ac (T5 - T4) + Ac (T5 - T2 ) + • • • + As (T5 - T5) = P5 , (19) 

where As is the thermal coupling parameter to the cold 
plate. We have programmed an algorithm that creates 
the thermal coupling matrix TCM. The thermal cou- 
pling equations can be written as 

TCMT = P + (20) 

We solve for the temperature distribution 



f = (TCM)"^ (p + As n 1) 



(21) 
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FIG. 12. a) diamond shaped resistor network b) mapping of 
the resistor network to a vector and a 2D matrix c) thermal 
cross couphng to nearest neighbors. 



We test for convergence by comparing the new temper- 
ature distribution with the old one. If the solution is 
within the convergence criteria, we increment the input 
current /in and calculate the next bias point, otherwise 
we recalculate the conductance vector G with the new 
temperature distribution and iterate the program. 
The RT curve is simulated in a similar manner. In this 
case the input current is fixed at a small value and the 
base temperature is changed by a small increment. 
The thermal cross coupling parameter Ac = 5.0 • 10 
was determined by considering the heat flow in the 
V2O3 film, only. The coupling parameter to the cold 
plate As = 4.25 • 10~^ ^ was obtained by optimizing the 
simulation of the IV curves with respect to the measured 
IV. The thermal conductance of a sapphire block with 
a ground area of 0.25 /im^ and a hight of 0.5 mm is 
approximately 7.5 • 10~^ ^. Here, it was assumed that 
sapphire at a temperature T = 150 K has a thermal 



7W 
K 



conductivity of 150 ^ [34. If one takes into account 
that edge cooling effects were neglected, the thermal 
coupling parameter As is in reasonable agreement with 
this estimation. 

In addition to the temperature distributions presented 
in the paper the voltage and current distribution was 
simulated as well. The current distributions depicted in 
Fig. ([8| and (10) are very similar to the temperature 
distributions presented in the main paper. 
The voltage distributions are shown in Fig. ([9| and ( 11 ). 
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